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The irregular polarity reversals of the Earth's magnetic field have attracted much interest during 
the last decades. Despite the fact that recent numerical simulations of the geodynamo have shown 
nice polarity transitions, the very reason and the basic mechanism of reversals are far from being 
understood. Using a paradigmatic mean-field dynamo model with a spherically symmetric helical 
turbulence parameter a we attribute the essential features of reversals to the magnetic field dynamics 
in the vicinity of an exceptional point of the spectrum of the non-selfadjoint dynamo operator. At 
such exceptional (branch) points of square root type two real eigenvalues coalesce and continue as 
a complex conjugated pair of eigenvalues. Special focus is laid on the comparison of numerically 
computed time series with paleomagnetic observations. It is shown that the considered dynamo 
model with high supercriticality can explain the observed time scale and the asymmetric shape of 
reversals with a slow decay and a fast field recovery. 

PACS numbers: 47.65.+a, 91.25.-r 

Introduction. 

There is ample paleomagnetic evidence that the axial dipole component of the Earth's magnetic field has reversed its 
polarity many times. The last reversal, the Brunhes-Matuyama transition, occurred approximately 780000 years ago. 
The mean rate of reversals varies from nearly zero during the Kiaman and Cretaceous superchrons to approximately 
5 per Myr in the present. Some observations suggest a pronounced asymmetry of reversals with the decay of the 
dipole of a given polarity being much slower than the following recreation of the dipole with opposite polarity 0, . 
Observational data also indicate a possible correlation of the field intensity with the interval between subsequent 
reversals [2, |3|- A third hypothesis concerns the bimodal distribution of the Earth's virtual axial dipole moment 
(VADM) with two peaks at approximately 4 x 10 22 Am 2 and at twice that value 0,H,E|- 

Although these reversal features are still controversially discussed in the literature, it is worthwhile to ask if and 
how they could be represented within geodynamo theory. With view on the recent dynamo experiments in Riga and 
Karlsruhe [7J one could also ask about the most essential ingredient for a dynamo experiment to exhibit irregular 
reversals in a similar way as the geodynamo does. 

In a recent paper we had shown that a simple mean-field dynamo model with a spherically symmetric helical 
turbulence parameter a can exhibit all three mentioned features of reversals. Interestingly, all of them are attributable 
to the magnetic field dynamics in the vicinity of an exceptional point [j| of the spectrum of the non-selfadjoint dynamo 
operator where two real eigenvalues coalesce and continue as a pair of complex conjugate eigenvalues. Typically, this 
exceptional point is associated with a nearby local maximum of the growth rate dependence on the magnetic Reynolds 
number. It is the negative slope of this curve between the local maximum and the exceptional point that makes even 
stationary dynamos vulnerable to some prevailing noise. This way, the system can leave the stable state and run 
towards the exceptional point and beyond into the oscillatory branch where the polarity transition occurs. 

A weakness of this reversal model was the apparent necessity to fine-tune the magnetic Reynolds number and /or 
the radial profile a(r) in order to adjust the operator spectrum in an appropriate way. In a follow-up paper |l0j it 
was shown, however, that this fine-tuning is not necessary in the case of higher supercriticality of the dynamo. It 
turned out that for increasing magnetic Reynolds number there is a strong tendency for the exceptional point and 
the associated local maximum to move close to the zero growth rate line were the indicated reversal scenario can be 
actualized. Although exemplified by the simple spherically symmetric a 2 dynamo model, the main idea of this "self- 
tuning" mechanism of saturated dynamos into a reversal-prone state seems transferable to other dynamos. Hence, 
reversing dynamos might be much more typical than what could be expected from a purely kinematic perspective. 

In the present paper we go one step further and compare paleomagnetic data of five reversals from the last two 
million years with numerical time series resulting from our model. It will be shown that it is again the strong 
supercriticality of the considered dynamo models that may explain the typical time scales of the observed asymmetric 
reversals. 
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FIG. 1: Magnetic field evolution for D = 0. With higher values of C, the field amplitude increases and the oscillation frequency 
decreases. Note also the growing anharmonicity (saw-tooth shape) of the oscillations and the transition to a non-oscillatory 
regime (as seen for C = 7.24). 



I. THE MODEL. 

We consider a simple mean-field dynamo model of a 2 type with a supposed spherically symmetric, isotropic helical 
turbulence parameter a flll | . The induction equation for the magnetic field B reads 

B = V x (qB) + (^oct) _1 AB , (1) 

with magnetic permeability [Iq and electrical conductivity a. For the Earth's core we will assume the diffusion time 
T diff = Hq<jR 2 to be ~ 200 kyr. 

The divergence- free magnetic field B can be decomposed into a poloidal and a toroidal components, according to 
B = —V x (r x V5) — r x VT. The defining scalars S and T are expanded in spherical harmonics of degree / and 
order m with expansion coefficients s/. m (r, r) and ti <m (r, t). For the envisioned spherically symmetric and isotropic a 2 
dynamo problem, the induction equation decouples for each degree / and order m into the following pair of equations: 

d 2 , v 1(1 + 1) , v 

d^( rsi > J3 - si + a{r,r)ti, (2) 

9 ( 9 . s , , d . A 1(1 + 1). . , , 

~dr [fr( rt ^~ a ( r > T >fr( rSl > ) [U - a(r,r)si] . (3) 

Since these equations are independent of the order m, we have skipped m in the index of s and t. The boundary 
conditions are dsi/ 'dr\ r —\ + (I + l)sj(l) = = 0. In the following we consider only the dipole field with / = 1. 

We will focus on a particular radial profile aki n (r) — 1.916 C (1 — 6 r 2 + 5 r 4 ), which is characterized by a sign 
change along the radius (the factor 1.916 results from normalizing the radial average of |a(r)| to the corresponding 
value for constant a). The motivation for this choice is that quite similar a(r) profiles had been shown to exhibit 
oscillatory behaviour ^1 ■ 

Saturation is ensured by assuming the kinematic profile akin (r) to be algebraically quenched by the magnetic field 
energy averaged over the angles which can be expressed in terms of Si(r, r) and t\(r, r). Note that this averaging over 
the angles represents a severe simplification, since in reality (even for an assumed spherically symmetric kinematic a) 
the energy dependent quenching would result in a breaking of the spherical symmetry 

In addition to this quenching, the a(r) profiles are perturbed by "blobs" of noise which are considered constant 
within a correlation time r corr . Physically, such a noise term can be understood as a consequence of changing boundary 
conditions for the flow in the outer core, but also as a substitute for the omitted influence of higher multipole modes 
on the dominant axial dipole mode. 



dsi 
dr 
dh 
dr 



3 




(a) ° 



(b) 



o 4.0 
o 3.0 
g- 2.0 
2 1.0 
0.0 



(c) 



2 3 4 
Time 





2 3 4 
Time 



0.4 
0.2 
; 0.0 
-0.2 
-0.4 







8 . 










' 6 5 4 
























" 3 1,2 • - 



(d) 



0.2 0.4 0.6 0.8 

r 




FIG. 2: Explanation of the field dynamics for C = 7.237 and D = 0. (a) Half an anharmonic oscillation with nine selected 
instants n, i = 1...9. (b) Instantaneous growth rates resulting from the instantaneous a(r,n) profiles, (c) Instantaneous 
frequencies, (d) Instantaneous poloidal field s(r,Ti). (e) Profiles a(r,Ti), compared to the unquenched a(r) (K). Note that the 
deformation of a(r) during a reversal is not very strong, (f) Details of (e) close to the origin. At the instants 6 and 7, the a(r) 
profile comes close to the unquenched a(r) (K). 



In summary, the a(r,r) profile entering Eqs. (2) and (3) is written as 

1.916(l-6r 2 +5r 4 ) 



a(r, t) 



C- 



i , p-l 2s\(rr) 1 / d{r Sl {r,r)) \ ,. 2 (\ 

Hi (r) + 6 (r) r 2 + 6 (r) r 3 + £ 4 (r) r 4 , 



(4) 



where the noise correlation is given by (^(t)^(t + n)) = D 2 (l — |ri|/r corr )6(l — \T\\/T corr )5ij. C is a normalized 
dynamo number, D is the noise strength, and E mag , n is a constant measuring the mean magnetic field energy. 



II. NUMERICAL RESULTS. 



We start with the noise-free case. Figure 1 shows the magnetic field evolution according to the equation system (2-4) 
for D = and different dynamo numbers C. The critical dynamo number is 6.78. The nearly harmonic oscillation 
for C = 6.8 becomes more and more saw-tooth shaped for increasing C, with a pronounced asymmetry between the 
slow field decay and the fast field recreation during the reversal. At C = 7.24 a transition to a steady dynamo has 
occurred. 

In order to understand this behaviour, we examine in Fig. 2 the evolution of the magnetic field within approximately 
half a period of the anharmonic oscillation for the particular value C — 7.237 at which the saw-tooth shape is already 
very pronounced. Figure 2a shows the time dependence of s(r = 1) during this half period, with the typical slow 
decay and the fast recreation of the field. This behaviour is analyzed in detail at 9 selected instants n (i = 1...9) for 
which the instantaneous fields s(r, r,) (Fig. 2d) and the corresponding a(r,Ti) (Fig. 2e, 2f) are depicted. The latter 
two plots reveal that a(r) undergoes only slight changes during the oscillation and that it comes very close to the 
unquenched, kinematic akin {r) (denoted by K) when the magnetic field is small in the middle of the reversal. 

It is instructive to plug the instantaneous a(r, Tj) profiles into an eigenvalue solver (for which the time derivatives 
on the left hand sides of Eqs. (2) and (3) are replaced by As(r) and Ai(r), respectively). Figure 2b and Fig. 2c show 



4 



o 
O 



K1 






I/O 

K2 












2 






3 


Sy ■ ' 




4 


' • V.. 




5 






6 vs? 

7 " • 

/ - 


Ej 


El, 2" 








■ i ><^* 


















/ ^6,7 , 







(a) 



4.5 



5.5 



6.5 



7.5 



7.237 




FIG. 3: The role of exceptional points for the understanding of reversals, (a) Instantaneous growth rate curves for the kinematic 
profile and for the quenched a(r,n) profiles. The Ei indicate, for each of the considered a profiles, the exceptional point, (b) 
Growth rates for the profiles a(r) = 1.916 C* C (1 - 6 r 2 + 5 r 4 )/(l + E mag (r)/E mag ,o) with C =6.78, 7.2, 10, 20, 50, 100, 
200. The label "7.2 (m)" refers to the maximally quenched a during the reversal. Note the remarkable move of the exceptional 
point well above the zero line and back to it (indicated by the thick dashed line). 



the resulting instantaneous growth rates and frequencies during the half-oscillation. Evidently, the reversal starts 
with a very slow field decay (slightly negative growth rate) which, however, accelerates itself and drives the system 
into an oscillatory behaviour (the frequency becomes different from zero for 2.0 < r < 4.5). Around t = 4, when the 
quenching of a is weak, the instantaneous growth rate reaches rather high values. Later we will see that for strongly 
supercritical dynamos these high growth rates result in a field dynamics that is much faster than what would be 
expected from the diffusion time scale. 

In Figure 3a we relate this magnetic field evolution to the existence of exceptional points of the spectrum. It shows 
the instantaneous growth rate which results from the individual a(r, n) profiles, and from the unquenched (kinematic) 
a(r) profile K. In addition to the growth rates at the actual dynamo number C = 7.237 (indicated by the dashed 
vertical line), we have plotted the whole growth rate curves in the interval between 4.5 < C < 7.5. At the exceptional 
point Efc of the kinematic aki n (r), the first eigenvalue branch Kl coalesces with the second branch K2 and both 
continue as a pair of complex conjugate eigenvalues. For all other curves, only the exceptional point is indicated by 
Ei whereas the branch of the second eigenvalue has been omitted. In this framework, a reversal can be described 
as follows: At the instant 1, the growth rate is located close to the maximum of the non-oscillatory branch which is 
slightly below zero. The resulting slow field decay accelerates itself, because the system moves down (instant 2) from 
the maximum of the real branch towards the exceptional point. Then the system enters the oscillatory branch (3-7). 
Finally the system moves back again (8,9) but with opposite polarity. 

The transition point between oscillatory and steady dynamos (C = 7.239) is characterized by the fact that the 
maximum of the non-oscillatory branch crosses the zero growth rate line. Beyond this point, the field is growing 
rather than decaying, leading to a stable fixed point somewhere to the left of the maximum of the non-oscillatory 
branch, and hence to a non-oscillatory dynamo (cf. the case C = 7.24 in Fig. 1). If noise comes into play it will soften 
the sharp border between oscillatory and steady dynamos. This means, in particular, that even above the transition 
point, the noise can trigger a jump to the right of the maximum from were the described reversal process can start 
(Fig. 2). _ 

Having identified the exceptional point and the nearby local maximum as the spectral features that are responsible 
for both the slow decay before and the fast recreation of the field after the polarity transition, one may ask now why 
the operator of the actual geodynamo should possess just such a special spectrum. 

A preliminary answer to this question was given in [lfj and is illustrated in Fig. 3b. Roughly speaking, highly 
supercritical dynamos seem to have a tendency to saturate in a state for which the exceptional point and its associated 
local maximum lie close to the zero growth rate line. Interestingly, this happens independently on whether the 
exceptional point in the original kinematic case was above the zero growth rate line or below it. Our example with 
otkinir) = 1-916 C (1 — 6 r 2 + 5 r 4 ) belongs to the second type. Starting from the kinematic akin{r) 1 for which the 
exceptional point is well below the zero line, it rises rapidly above zero to a maximum value, but for even higher C it 
moves back in direction to the zero line. 
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FIG. 4: Time series for some selected values of C and D. 



Without noise and for C > 7.239, the position of the local maximum above the zero growth rate line leads to a 
steady dynamo. The presence of noise will sometimes bring the actual growth rate below zero, and then the indicated 
reversal process can start. A number of time series for different values of C and D is depicted in Fig. 4. 

Figure 5 shows reversal details for four particular choices of C and D. For later comparison with paleomagnetic 
data we have exhibited five typical reversals, and their averages, from 80 kyr before the polarity transition until 20 
kyr after (except the case C = 8, D = 1 for which the reversal takes much longer). 



III. COMPARISON WITH PALEOMAGNETIC DATA. 



In this section we will validate if our model can be fitted to real paleomagnetic data. For this purpose we have 
used recently published material about five reversals which occurred during the last two million years [2|. Actually, 
the data shown in Fig. 6a have been extracted from Fig. 4 of 0. In all five curves, as well as in their averages, the 
asymmetry of the reversal process is clearly visible. The dominant features are a field decay over a period of 50-80 kyr 
and a rather sharp field recreation within 5-10 kyr. It has been an old-standing puzzle to explain in particular this 
fast recreation when a diffusive timescale of 200 kyr has to be taken into account. A possible solution of this problem 
is to assume a turbulent resistivity which is much larger than the molecular resistivity (e.g., by a factor 15-20 in \13§). 

The comparison of the real data with the numerically time series in Fig. 6b shows that this assumption of turbulent 
resistivity is by no means necessary. Apart from the slightly supercritical case C = 8, D = 1 which exhibits a much 
to slow magnetic field evolution, the other examples with C =20, 50, 100 and D = 6 show very realistic time series 
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FIG. 5: Field intensity (virtual axial dipole moment) variations across five selected reversals, and their averages, computed for 
different values of C and D. The time interval includes 80 kyr before (negative values) and 20 kyr after (positive values) a 
polarity transition. Note the quite different time scales for C = 8, D — 1 on one side, and for the other three examples on the 
other side. 



with the typical slow decay and fast recreation. As noted above, the fast recreation results from the fact that in a 
small intervall during the transition the dynamo operates with an nearly unquenched a(r) profiles which yields, in 
case that the dynamo is strongly supercritical, rather high growth rates. 



IV. CONCLUSIONS 

Based on former results in and |10| . we have shown that a simple but strongly supercritical a 2 dynamo model 
exhibits a number of features which are typical for Earth's magnetic field reversals, in particular an asymmetric shape 
and correct timescales for the field decay and field recreation. 

As it does not include the necessary North-South asymmetry of a we cannot claim that our model is an appropriate 
model of the geodynamo. However, recent papers by Giesecke et al. have shown that a) even in such realistic models a 
may exhibit a sign change along the radius 0] and b) that such models can also give rise to reversals . Therefore, 
it seems worthwhile to identify the indicated reversal scenario in this and in more complicated geodynamo models. 
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FIG. 6: Comparison of paleomagnetic reversal data and numerically simulated ones, (a) Virtual axial dipole moment (VADM) 
during the 80 kyr preceeding and the 20 kyr following a polarity transition for five reversals from the last 2 million years (data 
extracted from 0), and their average, (b) Comparison of the average curve of (a) with the four average curves of Fig. 5. 
The field scale for the numerical values has been fixed in such a way that the intensity in the non-reversing periods matches 
approximately the observed values. 
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